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We present a method for the measurement of a temperature differential across a single quantum 
dot that has transmission resonances that are separated in energy by much more than the thermal 
energy. We determine numerically that the method is accurate to within a few percent across a wide 
range of parameters. The proposed method measures the temperature of the electrons that enter 
the quantum dot and will be useful in experiments that aim to test theory which predicts quantum 
dots are highly efficient thermoelectrics. 



In the ongoing development of effective thermoelectric 
materials and devices, low-dimensional systems are of 
particularly great interest, because to optimize the per- 
formance of a thermoelectric, it is crucial to control the 
energy spectrum of mobile electrons [l], 0, d, 0, [BJ . De- 
vices for high-efficiency thermal-to-electric power conver- 
sion based on quantum dots defined by double barriers 
embedded in nanowires have been proposed [6|. Such 
systems have great advantages, because they select the 
energies at which electrons are transmitted 0, Q, and 
because nanowires can be contacted in highly ordered 
arrays @ with the potential for large-scale parallel oper- 
ation. 



In order to measure quantitatively the dependence 
of thermopower and energy-conversion efficiency on the 
transmission spectrum of a quantum dot, it is necessary 
to apply and determine accurately a temperature differ- 
ential across the dot. Traditionally for the thermoelectric 
characterization of mesoscopic devices such as quantum 
point contacts 10 1, quantum dots in 2DEG's UM . car- 
bon nanotubes 12 , [Tjj, [l4[ , and nanowires [ID, Il6| . an ac 
heating current generates a temperature differential that 
is measured in separate calibration experiments. Here we 
propose a technique that measures the actual electronic 
temperature differential across a quantum dot and does 
not require separate calibration. The basic concept is as 
follows: the change in current across a quantum dot in 
response to an applied heating voltage, Vkj is measured. 
This signal contains information about the electron tem- 
peratures at the source and drain, but it also depends on 
the dot's energy-dependent transmission function, r(e). 
However, one can obtain the necessary information about 
t (e) from conductance measurements. Together, these 
two measurements allow one to determine the source and 
drain temperatures separately. 



The two-terminal current through a quantum dot can 



be written [17|, [l 



[f s (e)-f d (e)]r(e)ds, 



(1) 



where — e^ sd + 1 are the Fermi-Dirac distributions 
in the nanowire's source and drain leads, respectively, 
and their arguments are £ Sj d = (e — fi s .d T e ^/2) /AsbT^cI. 
We assume the bias voltage, V, is applied symmetrically 
across the dot. For the case of a quantum dot or single- 
electron transistor (SET) with well-separated transmis- 
sion maxima as a function of gate voltage, Eq. (fl]) 
predicts the characteristic Coulomb blockade diamonds 
which appear in the differential conductance, G, as a 
function of bias voltage and gate voltage. 
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FIG. 1: The heating setup and the temperature landscape. 
The source contact is warmed with a voltage-balanced heating 
current and can be biased for thermocurrent measurements. 
Electron transport through the quantum dot is determined 
by the local temperatures of the source and drain sides of the 
dot, T s , d = AT s , d + T . (Color online) 

In a typical experiment, an ac heating current is used 
to modulate the temperature Tjr of an ohmic contact 
at one end of a nanowire (taken here to be the source 
contact, Fig. [5)19] with amplitude ATh with respect 
to the unperturbed device temperature, To. We are in- 
terested in the associated electronic temperature rises, 
AT Sj d = T s ,d — To, in the immediate vicinity of the quan- 
tum dot (see Fig. []}. In the case of strong electron- 
phonon interaction (for example near room temperature) 
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the electronic temperature will drop linearly along the 
nanowire, and AT S = ATd if the quantum dot is short 
compared to the wire. At low temperatures, however, 
where electron-phonon interaction in the nanowire is ex- 
pected to be weak, AT H > AT S > AT d > 0, and AT S 
and ATd need to be measured. 

Assuming an ac heating voltage, Vh = Vo cos(wt), the 
temperature rises on the source and drain sides of the dot 
can be written AT Sj d = /3 s ,d^H > wnere Ps,d are unknown 
constants and 7 can take on various values depending on 
the type and strength of electron-phonon interaction in 
the heating wire [20]. Here we assume a short heating 
wire and Joule heating, and therefore 7 = 2 [2(|. In 
this regime, by an application of the chain rule, the rms- 
amplitude of the ac temperature rises can be written 
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In an experiment, one can measure 01 /0 (Vjf) , the 
frequency-doubled response to the ac heating voltage. 
The differential thermocurrent, 
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cannot be measured directly. However, we will show 
that it can be obtained in good approximation from con- 
ductance measurements. 

Under bias conditions, where the source (drain) elec- 
trochemical potential is near a well-defined transmission 
resonance of the quantum dot, while the drain (source) 
is several k^T away from the next resonance, the second 
derivative of the current is 
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A key observation is that the integrands in Eq. and 
Eq. ( [5]) are qualitatively very similar: 
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This approximation holds for all £ Si d, because 2/ Si d — 1 
limits to — £ s ,d/2 when £ Sj d is small, and df s ^/d£, s ^ goes 
to zero in all other cases. With this approximation, we 
can combine Eqs. ([3]) and : 
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where A is a unitless scaling factor introduced during in- 
tegration. In this way, all the information about r (e) 
needed to determine 01 /dT s ^ is accounted for by mea- 
suring 2 I/0V 2 . 

Substituting Eq. ([6|) into Eq. ([2|) and solving for AT Si d 
yields our final result, 
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FIG. 2: A plot of dI/dT H (red) and (2fc B /e) 2 AT H d 2 I/dV 2 
(blue) as a function of bias voltage at a fixed gate voltage, as 
indicated by the horizontal dashed line in Fig. [3] The temper- 
ature rise (green) is calculated via Eq. @. Insets: The posi- 
tion of the resonant tunneling energy of the dot with respect 
to the electrochemical potentials in the leads. In this model, 
we used a transmission function consisting of Lorentzians with 
a FWHM of T = 0.5 meV equally spaced by 5 meV and 
AT H = 300 mK, AT = 0.8AT H = 240 mK, AT d = 0.2AT H 
= 60 mK, and T = 230 mK. (Color online) 
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FIG. 3: The calculated temperature rise as a function of bias 
and gate voltages for a single Coulomb blockade diamond. In 
the regions indicated, either the source or drain electrochemi- 
cal potential, but not both, is within a few A?bT of a resonant 
energy of the dot. Our assumptions are fulfilled in these re- 
gions, and the simulation produces temperature plateaus pre- 
dicting the correct AT s ,d- Fig. [2] shows a slice at 36.75 mV 
gate voltage, as indicated by the dashed line. (Color online) 



which shows that approximations of AT S and ATd can be 
obtained from measurement of 2 I/0V 2 and 0I/0(V^) 
and knowledge of T). 

To illustrate the qualitative similarity of Eqs. ([3]) and 
(HI), we show numerical calculations of the two in Fig. [5] 
taken at the gate voltage indicated by the horizontal 
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FIG. 4: The calculated temperature rises as a function of 
ATh . The calculated values agree well with the expected 
values (solid lines) up to nearly 10 To- Inset: The percent 
error as a function of T, the full width at half max (FWHM) 
of the transmission function, r (s) in Eq. fl} . The error is 
within 1% over an order of magnitude in F. Here AT S is 240 
mK and ATd is 60 mK. (Color online) 

dashed line in Fig.[3]at 36.75 mV. The left inset of Fig. [2] 
illustrates that, in this example, when the bias voltage 
is negative, the source temperature is the only temper- 
ature affecting the current through the dot; therefore, 
dl 7<9Th = dI/dT s in this bias configuration. In the op- 
posite configuration, dl /9Th = dl/dT^, as shown in the 
right inset of Fig. O Fig. [3] shows Eq. Q as calculated 
from modeled data of Eqs. ([3]) and (|4|) across an en- 
tire Coulomb blockade diamond (indicated by a white, 



dashed line), and a slice through that diamond is shown 
as green symbols in Fig. [2l In regions along the diamond 
ridges (circled areas in Fig. [3]) — where one, and only one, 
of the two electrochemical potentials in the source or 
drain is within a few k^T of a transmission resonance — 
Eq. (O yields consistent values in accordance with the 
assumptions that allow us to write Eq. ((4|. In all other 
regions, Eq. ((¥]) is not valid, not even approximately, be- 
cause it only accounts for one Fermi-Dirac distribution. 

The use of this method requires knowledge of the ap- 
propriate scaling factor A, defined in Eq. [|5J which needs 
to be determined numerically. For the particular mod- 
eling parameters used here (see caption of Fig. [2]), we 
found A = 0.304 by averaging Eq. ([7|) over the voltage 
range from the peak value of dl/dT^ to 20% of its peak 
value, where the signal-to-noise ratio in an experiment 
should be largest. Note that for different parameters, 
A will differ, but it is insensitive to typical experimen- 
tal variations in T and ATh- For example, in Fig. 2J 
we show that the use of the same A = 0.304 (calculated 
for r = 0.5 meV) yields errors in AT S and ATd of only 
1% when r is varied over nearly an order of magnitude 
around T = 0.5 meV (inset of Fig. [4]) and only a few per- 
cent for ATh up to almost 10 To- To put this small error 
into context, note that the local temperatures T s and Td 
can be defined only over a distance of about an inelastic 
scattering length, such that an accuracy of less than a 
few percent is not necessarily physically meaningful. 
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